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Abstract: In living cell or its nucleus, the motions of molecules are complicated due to the large crowding and expected heterogeneity of the 
intracellular environment. Randomness in cellular systems can be either spatial (anomalous) or temporal (heterogeneous). In order to separate both 
processes, we introduce anomalous random walks on fractals that represented crowded environments. We report the use of numerical simulation 
and experimental data of single-molecule detection by fluorescence fluctuation microscopy for detecting resolution limits of different mobile 
fractions in crowded environment of living cells. We simulate the time scale behavior of diffusion times Td(t) for one component, e.g. the fast 
mobile fraction, and a second component, e.g. the slow mobile fraction. The less the anomalous exponent a the higher the geometric crowding of 
the underlying structure of motion that is quantified by the ratio of the Hausdorff dimension and the walk exponent d f/d w an d specific for the 
type of crowding generator used. The simulated diffusion time decreases for smaller values of a ^ 1 but increases for a larger time scale T at a 
given value of a ^ 1. The effect of translational anomalous motion is substantially greater if a differs much from 1. An a value close to 1 
contributes little to the time dependence of subdiffusive motions. Thus, quantitative determination of molecular weights from measured diffusion 
times and apparent diffusion coefficients, respectively, in temporal auto- and crosscorrelation analyses and from time-dependent fluorescence 
imaging data are difficult to interpret and biased in crowded environments of living cells and their cellular compartments; anomalous dynamics on 
different time scales 1 must be coupled with the quantitative analysis of how experimental parameters change with predictions from simulated 
subdiffusive dynamics of molecular motions and mechanistic models. We first demonstrate that the crowding exponent a also determines the 
resolution of differences in diffusion times between two components in addition to photophyscial parameters well-known for normal motion in 
dilute solution. The resolution limit between two different kinds of single molecule species is also analyzed under translational anomalous motion 
with broken ergodicity. We apply our theoretical predictions of diffusion times and lower limits for the time resolution of two components to 
fluorescence images in human prostate cancer cells transfected with GFP-Ago2 and GFP-Agol. In order to mimic heterogeneous behavior in 
crowded environments of living cells, we need to introduce so-called continuous time random walks (CTRW). CTRWs were originally performed 
on regular lattice. This purely stochastic molecule behavior leads to subdiffusive motion with broken ergodicity in our simulations. For the first 
time, we are able to quantitatively differentiate between anomalous motion without broken ergodicity and anomalous motion with broken 
ergodicity in time-dependent fluorescence microscopy data sets of living cells. Since the experimental conditions to measure a selfsame molecule 
over an extended period of time, at which biology is taken place, in living cells or even in dilute solution are very restrictive, we need to perform 
the time average over a subpopulation of different single molecules of the same kind. For time averages over subpopulations of single molecules, 
the temporal auto- and crosscorrelation functions are first found. Knowing the crowding parameter a for the cell type and cellular compartment 
type, respectively, the heterogeneous parameter y can be obtained from the measurements in the presence of the interacting reaction partner, e.g. 
ligand, with the same a value. The product a ' y = y is not a simple fitting parameter in the temporal auto- and two-color crosscorrelation 
functions because it is related to the proper physical models of anomalous (spatial) and heterogeneous (temporal) randomness in cellular systems. 
We have already derived an analytical solution for J in the special case of y = 3/2 . In the case of two-color crosscorrelation or/and two-color 
fluorescence imaging (co-localization experiments), the second component is also a two-color species gr, for example a different molecular 
complex with an additional ligand. Here, we first show that plausible biological mechanisms from FCS/ FCCS and fluorescence imaging in living 
cells are highly questionable without proper quantitative physical models of subdiffusive motion and temporal randomness. At best, such 
quantitative FCS/ FCCS and fluorescence imaging data are difficult to interpret under crowding and heterogeneous conditions. It is challenging to 
translate proper physical models of anomalous (spatial) and heterogeneous (temporal) randomness in living cells and their cellular compartments 
like the nucleus into biological models of the cell biological process under study testable by single-molecule approaches. Otherwise, quantitative 
FCS/FCCS and fluorescence imaging measurements in living cells are not well described and cannot be interpreted in a meaningful way. 
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Walks (CTRW) on fractal supports, cellular crowding, Cytoplasmic Assembly of Nuclear RISC, ergodicity, FCS, FCCS, 
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1. INTRODUCTION 

In living cells, the abundance of biomacromolecules in 
the cytoplasm and nucleus causes a state of molecular 
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crowding [1-4]. Translational motion in crowded envi- 
ronment is not well represented by normal, Brownian 
diffusion [5]. Translational diffusive measurements by fluo- 
rescence fluctuation microscopy (FFM) like fluorescence 
correlation spectroscopy (FCS) are difficult to interpret in a 
crowded environment of a living cell because conventional 
analysis of fluorescence fluctuation time series such as auto- 
and crosscorrelation function assumes that the spatial domain 
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available for diffusive motion is unlimited [6]. This assum- 
ption is invalid for a crowded environment of a biological 
cell. The overall translational motion of molecules in living 
cells is better described by distances between collisions that 
are not independent and can take different times to perform. 
The temporal autocorrelation and the temporal crosscorre- 
lation of the signal with itself reveals the time required for a 
molecule to diffusion through an open observation volume 

A V. Diffusive motion is the macroscopic result of random 
thermal motions that occur on a molecular level, and the 
collision between solute particles/molecules and solvent 
molecules experiencing thermal movement are responsible 
for this phenomenon. Along this line, translational anoma- 
lous diffusion is interesting to explore because it is a 
prototype of subdiffusion in biological cells and their 
compartments like membranes [7]. A common way to study 
deviations from normal Gaussian diffusion is to plot the 
mean square displacement (MSD) of the molecules and 
particles, respectively, in ^-dimensional space. In systems, 
where molecules move freely via 3D Brownian dynamics, 
the MSD is proportional to time and given by 

r m\ = 6Dt- D and t are the diffusion coefficient and 



time. Any deviation from linearity corresponds to 
translational anomalous diffusion [8] with 



r (t) 

and a ^ 1. We focus on subdiffusion, i.e. 0 < a < 1. While 
the MSD is used to classify a process as subdiffusion, it does 
not provide any information on the physical and biological 
mechanism underlying the subdiffusion. A physical correct 
propagator for subdiffusive temporal evolution based on the 

instantaneous diffusion coefficient D imt (t) = —.Y-t a ~ l was 

6 

first described, producing the correct power law scaling of 
MSD versus time and rigorously solved the extended 
diffusion equation from Fick's law and the continuity 
equation [9]. With that propagator, the correct FCS 
autocorrelation function for translational anomalous 
diffusion was derived [9], where the temporal 
autocorrelation defines the time scale X of diffusion with the 
characteristic diffusion time X D across the 3D observation 
volume AV for »2-photon excitation and ^-dimensional 
motion as 



Tn = 



f > 


l/a f 






4 ■ m ■ r/2 • n 




V J 


V 



or 



a-(o,_ 



(1) 



m is the radial waist of excitation laser and n the 

*"x-y app 

measured or apparent diffusion coefficient, e.g. in FCS. With 
this clarification, it becomes more feasible to unambiguously 
report mobility data in terms of either a well-defined diff- 
usion time T D or as time-dependent diffusion coefficients [9] . 
Two-color crosscorrelation analysis (FCCS) for one-photon 
excitation only requires in Eqn. (1) that 
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where CO and CO are the radial waists of the green 

x-y,g w x-y,r & 

(or blue) and red excitation lasers for the measurement of the 
two-color labeled molecular species gr. The most commonly 
published form of the time-dependent 3D diffusion 
coefficient is defined as [7] 

A P p(0 = --' 
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such that the MSD can then be written as 



r (x)) = 6-D app (r)-r. 



(4) 



The apparent 3D diffusion coefficient D app (T) is found 

at time T from the slope of a straight line between the origin 
and MSD(t) in the linear plot of MSD(t) versus T [9]. Using 
FFM and FCS to investigate translational anomalous 
mobility in living cells requires a through understanding of 
the theoretical basis of fitting models and the physical 
significance of associated fitting parameters as we shall 
demonstrate in this original article. 

In FFM, movement of fluorescent molecules through 
subfemtoliter observation volume AV generates temporal 
fluctuations in the fluorescence emitted across the 3D 
observation volume. We demonstrated by simulations and 
experiments that at the single-molecule level there is only 
one molecule at a time inside AV or no molecule in dilute 
solution-phase [10]. With the average molecule number N = 
0.0055 for AV = 0.21 fL measured with 24-nm nanospheres 
at 635 nm laser excitation and N = 0.0052 for AV = 0.14 fL 
measured with 100-nm nanospheres at 470 nm laser excita- 
tion in the real fluorescence fluctuation time series, we 
reached the lower single-molecule detection limit for dilute 
solutions. The contribution of two simultaneously fluore- 
scing molecules to the detected signal during the measure- 
ment vanished for N « 1 [10], and the single-molecule 
detection regime came close to the average N = 0.048 and N 
= 0.0057 molecules and particles, respectively, per AV [11]. 
In this paper, we analyze translational anomalous, subdiffu- 
sive motion without broken ergodicity and with broken ergo- 
dicity in the 3D-measurement set that includes an ensemble 
average of single molecules, or a single molecule, an obser- 
vation volume AV and a local crowded or/and heteroge- 
neous environment. The article aims to remove some wide- 
spread confusion regarding the interpretation of subdiffusive 
measurements in crowded and heterogeneous environment of 
living cells by time-dependent fluorescence fluctuation 
microscopy. 

We apply our theoretical predictions of diffusion times 
and lower limits for the time resolution of two components 
to fluorescence images in human prostate cancer cells (PC-3) 
transfected with fusion proteins of green fluorescence protein 
(GFP) and argonaute proteins (Agol, Ago2). Experiments in 
living cells and their compartments are challenged by the 
fact that diffusion is slower and heterogeneous compared to 
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in vitro systems of dilute solutions [6]. The cell response 
depends on a spatially organized system of biomacro- 
molecules like proteins, RNAs, DNA. To sustain cellular 
activity, the spatial crowding of biomacromolecules is main- 
tained within the intact cellular architecture [12-14]. We first 
address the problem how translational anomalous, subdif- 
fusive motion dominates the time-dependent behavior of an 
ensemble average of single molecules or a single molecule. 
We shall show that ignoring the observed effects of crow- 
ding and heterogeneity may produce results that cannot be 
interpreted in a meaningful way. 

2. METHODS 

2.1. Random Walks on Fractals 

There is a pronounced interest in clarifying the behavior 
of random walks (RW) on fractal supports [15]. The appro- 
aches that we use are random walks on a fractal support and 
continuous time random walks (CTRW) [8, 16]. Both appro- 
aches aim on the clarification of the interrelation of random 
processes connected with real diffusion phenomena observed 
in experiments. The majority of work found in literature, 
however, is restricted to one dimensional or two dimensional 
systems [8, 17]. Only a few articles deal with the real three 
dimensional problem [15]. In uniform Euclidean systems, the 

mean-square displacement ( r ) of a random walker is 



proportional to the time t for any number of spatial 
dimensions d . However, in disordered systems, this law is 
not valid in general, and the diffusive law becomes 
anomalous [18, 19] 



= r -t a 

1 a 1 ' 



(5) 



which in addition with the first moment, the mean, is 
discussed as a finger print of normal (a = 1) and anomalous 
(0 < ce < 1 ) diffusion [17]. It is well known that diffusive 
motion on a fractal support is subdiffusive with the MSD in 
time t given by Eqn. (5). The exponent a is related to the 
fractal dimension of the walker d =21 a [20]. The slowing 

down of the motion is caused by the delay of the diffusing 
molecules in the dangling ends, bottlenecks and backbends 
existing in the disordered structure. On the other hand, Eqn. 

/- 2 \ 

(5) is based on the existence of the term ( r ) , which is for 



an ergodic system either the ensemble mean over all random 
trajectories or the time average of a single trajectory [8]. 
Thus, the first approach in this article is to examine the 
motion on a fractal support as an ensemble of tracks carried 
out on it. Two and three dimensional examples of tracks of 
random walk on fractal support are shown in Fig. (1). In a 
simple discrete random walk, the walker advances one step 
in unit time. Each step is taken to a nearest neighbor of the 
site. The simple random walker can step with equal 
probability to any of the nearest neighbor sites that belongs 



to the grid. There are several ways of assigning transition 
probabilities for stepping from site to site [15]. The motion 
of the Brownian molecule is generated with a random 
number generator delivering pseudo-random numbers used 
for the steps in all three spatial directions. We generate a 
random Brownian walk by randomly selecting steps in the 
three coordinate directions. The three coordinate directions 

are generated by a permutation of the vector v = (0,0,1) so 

that a set of orthogonal vectors S is generated. 
Mathematically this means we use the basic set of orthogonal 
unit vectors in a Cartesian coordinate system as the basis of 
our calculations 



S = \ i,j,k\ i,j,ke. 



A 



= 1, 



= 1, 



= 1 



(6) 



This set of permuted vectors is extended in all directions 
positive and negative by the following unification of basis 
sets 



S =Su-S= i,j,k,-i,-j,-k\. 



(7) 



Introducing the random function R k which selects the 
direction with equal probability randomly from our basis set 

S*, we create the Brownian track B [ ro,r ] by a sum of 



(8) 



independent vectors 

B track :=^+IX(S*) = B J ro,r\, 



A simple random walk is statistically self-similar [15, 8]. 
The random steps are distributed according to a probability 

function P ( r 1 • In the limit that n » 1 , P I r 1 tends to a 



Gaussian distribution. This is a simple result of the Central- 
Limit Theorem. A direct consequence of the self-similarity is 
that the simple random walk is a statistical fractal. Upon 

dilation of space by a factor of A 12 , the number of steps 
increases by a factor of A . Therefore, the fractal dimension 
of a random walk is d w = ln(A)/ ln^l 1/2 )= 2 • It is a curious 

result that random walks performed on disordered but 
statistically self-similar structures are still self-similar 
themselves, exactly as in Euclidean space. The quantitative 
difference is that the usual diffusion exponent, d =2 is no 

longer equal to 2 and diffusion becomes anomalous. In order 
to generate an anomalous random walk, we use the set S* 
and the random function R t restricted to a fractal structure 

which is produced on a virtual grid. Examples of such tracks 
are shown in Fig. (1) for two and three dimensions. This 
approach can be easily generalized to higher dimensions if 
the basis set S is generalized to the corresponding span of 
the space under consideration. How the fractal support is 
generated will be discussed below. 
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Fig. (1). Representation of fractal tracks in two and three dimensions. Left: a track on a Sierpinski gasket. Right: three different trajectories 
on a generalized 3D Sierpinski carpet. The left panel shows that the 2D random walk on fractal support is traversing the fractal geometry 
densely in some parts of the fractal. The right panel shows that a 3D random walk on the fractal support is examining the space on the fractal, 
however the fractal support is not obvious for these short tracks with a total number of steps n = 10 4 . 



The finger print of RW is its mean square displacement. 
The MSD on a 3D random walk on fractal support is 
compared with the 2D-RW in Fig. (2). The MSD in all cases 
is determined by the ensemble average of M tracks of 

length n using the random track B n ^ r o , r J 

If ergodicity is given an equivalent representation can be 
formulated for the MSD, a moving average over a single 
trajectory of length n 

(7j =^-J o T " 4 ^2(f+A)-l(t)j-^a+A)-B(o^, ( 10 ) 

where B(t) is the continuous representation of B n f ro.rj 

and A is the lag time for the continuous random walk of 
length T. T is the measurement time in the real experiment. 

2.2. Generation of Fractals 

The term fractal was introduced by Mandelbrot [21] and 
is nowadays used to distinguish irregular structures in nature 
that are self-similar on certain length scales [22]. Fractals are 
showing up almost everywhere, for examples as coastlines, 
cells, human lungs, or percolation clusters [23]. The men- 
tioned examples are so called statistical fractals that, in 
contrary to deterministic fractals, show the self-similar beha- 
vior only on a limited length scale [23-27]. However, 
deterministic fractals show a self-similar behavior on all 
length scales. One of the most prominent examples of a 
deterministic fractal is the Sierpinski gasket in 2D [8] and 
the Sierpinski carpet (Menger sponge) in 3D. To construct 
the Sierpinski gasket we start with a square which is divided 
into nine equal subsquares, where the central square is 
removed. This first step of the construction is called the 
generator of the 2D Sierpinski gasket. If we continue with 
the construction in the same way for each of the remaining 



subsquares the second iteration step of the fractal is 
generated [8]. Continuing to infinity, we will get a structure 
with a self-similar behavior on all scales containing many 
holes arranged in a symmetric way. The same procedure can 
be used to construct the 3D version of the Sierpinski gasket 
which is also self-similar and has the appearance of a 
sponge. However, the 2D and 3D construction are not 
restricted to this simple pattern. In fact, it is possible by the 
same method to generalize the Sierpinski gasket and the 
carpet to a different structure, if we not only delete one 
element in the generator but, instead, allow the deletion of 
more than one element. This, of course, results into a great 
variety of generalized Sierpinski patterns introducing a 
variation of the gasket and the sponge. In our examinations, 
we will restrict us to generalized Sierpinski carpets (GSC), 
which delete not more than half of the elements of the 
generator in d dimensions. 

We begin by considering a class of fractal subsets of R'' 
formed by the following generalization of the construction of 
the Cantor ternary set. Let d > 2 and let 

S o = [°> 1 f =[1,0]X[1,0]X...X[0,1] the d -dimensional Carte- 
sian product. Let v > 3 be an integer and divide S 0 into v d 

equal subcubes. Next remove a symmetric pattern with o 
elements of subcubes from S 0 and call what remains S l ■ 

Now repeat the procedure: divide each subcube that is 
contained in S l into v d equal parts, remove the same 
symmetric pattern from each as was done to obtain S l from 
S 0 , and call what remains S 2 ■ Continuing in this way we 
obtain a decreasing sequence of (closed) subsets of [0,1]'' • 

Let S = n S ; we call S a generalized Sierpinski carpet 

(GSC) or simply, a carpet. We may regard the set S as 
idealized models of a region with diffusion obstacles of 
many different sizes. The standard SC is the GSC for which 
d = 2, V = 3 , and S l consists of S 0 minus the central 
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Fig. (2). The mean-square displacement of a normal random walk (RW) on a 2D regular square grid (left panel) and the MSD of a 3D 

h 2 \ 

random walk on a generalized Sierpinski carpet (right panel). For large track lengths n, it is obvious that the linear behavior / y ) ~ n 

h 2 \ 

changes to the anomalous ( r ) ~ n" ■ The case shown with a = 0.98 was calculated for a generalized Sierpinski carpet. The number of 



tracks used to derive the MSD are in both cases M = 10 3 . In both plots the straight line represents the linear behavior. For anomalous 
diffusion, a significant deviation corresponding to the law n a is obvious. 



square. Let ji = v — o be the number of subcubes 
remaining in S l , and let 



d, = ln(//)/lri(v). 



(11) 



Hence, the Hausdorff dimension of 5 is df. An example of 
a GSC in R 3 , is depicted in Fig. (3), where d = 3 and v = 3 
with \i = 15. 




Fig. (3). Construction of a generalized Sierpinski carpet with V = 3 
and U = 15 (12 empty sites of deleted cubes). The right graph shows 
part of the 4th iteration of the fractal support. In the simulations of 
experimental data we actually used the 12th iteration of the fractal 
which is not shown here due to space requirements. 

This generation of a fractal is based on the 
renormalization of the whole structure and can be used 
efficiently to generate a fractal support on an infinite space. 
Since ji<v i for d = 3 is the number of occupied sites that 

can be arranged on V 3 sites in an independent way, a large 
number of arrangements is possible. In fact, if we count all 



the different arrangement we will find p = 

[n 

arrangements; e.g. \i = 22 and v = 3 allows p = 80730 
possible patterns. This means that the arrangement of the 
occupied cubes determine the various patterns, which can be 
created each having the same Hausdorff dimension d f ■ 

However, the partition of the space will be different for each 
of the unique generators. It becomes obvious due to the 
binomial behavior of arrangements that there exist a large 
variety of generators possessing all the same dimension d f ■ 

After discussing how fractals are generated on latices, we 
perform random walks on these lattices. As these fractals 
have holes on every length scale due to their construction 
procedure, the diffusion on such structures is slowed down 



= r-t a =r n 



(12) 



The walk exponent d is usually greater than two for 

fractals. The relation a = 2/d < 1 is valid for any Euclidean 

dimension. For a random walk on a fractal support, t is 
equivalent to the number of steps n . 

A random walk (RW) along the occupied sites gives the 
walk exponent d from Eqn. (12) [28]. The value of d is 

determined not only by d f but also by other details of the 
generator. The deviation between d w and d f are discussed 

using the term lacunarity to explain the differences [21, 28, 
29]. A relation between the Hausdorff dimension and the 
walk exponent using the equivalence of the vibration and 
diffusion problem was derived by Alexander and Orbach 
[18] 
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2 



(13) 



where d s is the spectral dimension. The great interest in the 

conjecture results not only from the universal value assigned 
to d = 4/3 , but also from the fact that it provides a 

relationship between dynamic and static exponents. Indeed, 
by assuming c/ s = 4/3 = constant in the exact relation, one 

immediately obtains d w = 2d f /d s , characterizing the 
dynamic process of diffusive motion, in terms of d f , which 
describes the static geometrical scaling of the mass of the 
molecule. While the exact determination of d s seems to be 

a hard problem, it is quite easy to obtain certain bounds, and 
we have in particular that 1 < J s < d f < d ■ 

The majority of random walk algorithms in use for RW 
on fractal supports have huge memory requirements, as a 
complete fractal with a fixed iteration depth is stored. In our 
algorithm, we do not store a generated pattern of the fractal 
up to some definite stage [8, 17]. We only supply the 
generator as input. The random walker is set on some site 
and it tests whether each site it arrives at is an allowed site, 
as it goes along. This kind of walk generation is known as 
the blind ant approach [30]. The actual procedure is as 
follows: The walk can start at any site of the underlying 
virtual lattice. To check whether a site is accessible, the first 
step is to identify the iteration stage the point belongs to. For 

an V 3 grid, a point having either an X (<) -coordinate 

k— 1 k j 

(i = 1,2,3) between V and V belongs to the kth 
iteration stage of the fractal. In the k th-stage coarse-grained 

k— 1 

pattern with units of size V , it is checked whether the 
block containing the site matches an accessible site on the 
given generator. If found accessible, the corresponding point 
in the next lower stage, i.e. (A: — 1), is ascertained. In this 
way, the point is successively scaled down until it reaches 
the first stage. In general, in the k th stage, the equivalent 
coordinates (rf ) are given by the integer parts of 

xf/v*" 1 with / = 1,2,3. If feVf.xf ) matches an 
allowed site, the coordinates carried over to the next stage 
are x f_, = mod(xf,v i " 1 ) ( = 1,2,3- If an equivalent 

coordinate of any stage does not match the list of accessible 
sites, the site under consideration is blocked, only those sites 
surviving up to stage 1 are accessible. If the point corres- 
ponds to a blocked site, at any stage of the process, it is inac- 
cessible. This procedure of coarse graining the grid corres- 
ponds to a renormalization of the lattice [15]. It is useful to 
know conditions under which scaling exponents a and 
values of T a are generated through the random walk on the 

3D Sierpinski carpet. As exemplified in Table 1, a and T a 

are a consequence of the chosen generator of the 3D 
Sierpinski carpet. 



2.3. 3D Measurement Set 

This original research article deals with a 3D 
measurement set which includes an ensemble of single 
molecules, or a molecule, an observation volume AV and a 
local environment [10]; in contrast to our previous paper 
[10], in this article the local environment is the crowded and 
heterogeneous environment like a biological cytoplasm or 
the cellular compartment of a nucleus and not a dilute 
solution. As explained above and depicted in Fig. 1, single- 
molecules show translational anomalous diffusive motion, 
which we simulate by fractal dynamic behavior of the 3D 
random Brownian movements. In fluorescence fluctuation 
microscopy, no information can be obtained about the single 
molecules as long as the molecule dwells outside the 3D 
observation volume AV. The dynamics of the molecule 
remains hidden unless it is the dynamics of the molecule 
itself that causes the change in the molecule number 
fluctuations across AV. For single-molecule fluctuation 
counting, we introduced the scalar function T\ (t) for the real 
experimental situation of measuring a single molecule at a 
time [10] 

fl inside AV, emitted fluorescence detected, 

7](0= ( 14 ) 
[0 outside AV, no florescence detected. 

Consequently, there is a simple in and out jump of the 
molecule or an on and off of the signal constrained by AV as 
shown in Figs. (2 and 4). This represents the real molecule 
number fluctuation for the measurements. 

3. RESULTS AND DISCUSSION 

A conceptual framework for crowded and heterogeneous 
environments like biological cells and their cellular 
compartments is translational anomalous, subdiffusive 
motions of biomacromolecules due to geometric and/or 
energetic disorder [8]. The motions of proteins, RNAs or 
DNA can be hindered either by molecular crowding modeled 
by spatial, geometric restraints of fractals or by energetic 
landscapes/chemical binding represented by continuous-time 
random walks with heavy tails. This has been coined a 
coexistence of an ergodic anomalous diffusion and a 
nonergodic process of dynamics of macromolecules in living 
cells and their complex environment. Deciphering the 
molecular nature of the coexistence of ergodic anomalous 
process and nonergodic behavior constitutes an important 
quest of the field. 

3.1. Translational Anomalous Motion with Unbroken 
Ergodicity: Bridging Single-Molecule Approaches with 
Ensemble Averages 

We first develop analytical models for the simulation of 
3D subdiffusive motion. Widely used experimental methods 
such as fluorescence fluctuation spectroscopy and imaging 
like fluorescence correlation spectroscopy reveal only 
apparent values of the diffusion coefficients as well as a 
values. As shown in Table 2, the developed fractal 
generators yield realistic scaling exponents a and T a 
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Fig. (4). The measurement system: 3D translational anomalous 
diffusive motion within the observation volume AV = 0.14 fL (in 
pink color). Simulation steps n=10000. 
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values of the 3D random walks on the their structures. 
Molecular crowding means that the possibilities of molecule 
motion are geometrically restricted by the local cellular 
environment where the molecule motion takes place. The 
less the exponent a the higher the geometric crowding of the 
underlying structure of motion, which can be quantified by 
the ratio of d f ld n . as depicted in Fig. (5). The ratio is 

specific for the type of generator used. 

For a confocal observation volume AV of0.14fL(l fL = 
10~ 15 L) with a radial waist (O x _ v =196 nm of the excitation 
laser and a half-length of co = 0.478 nm, we measured an 
average diffusion time T D , = 0.628 ms over all time scales T 
across this 3D observation volume AV for the fast mobile 



Anomalous Diffusive (Subdiffusive) 3D Motion on 3 Fractal Supports. Type of Fractal: 3D Sierpinski Carpet. The Fractal 
Generators are Depicted in the Model. The Simulated 3D Random Walks on Each Fractal Support are Shown in the Right 

Panels. Simulation Steps n = 10000 . Different a and Y a Values were Obtained from the Double-Logarithmic Plot of the 

Simulated MSD(f) versus t Generated from the Different Fractal Supports. For Spatial and Time Steps of Unity, the 

Following a and T were Simulated from Top to Bottom: a = { 0.863, 0.766, 0.770 } , V = { 0.158, 0.355, 0.332 } 
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Table 2. Simulation of Values for a, T and Time Steps At with the Indicated Fractal Generators by Fixing the Spatial Step 
Lengths of the 3D Random Walks Ax to 10 nm on the Crowded Structure 

Generator 
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Fig. (5). Relation between the fractal ratio d f ld as a function of 

CL . The fixed ratio of Alexander and Orbach (AO) with 
dJd w = 3/2 is also shown as a dashed line (top of the graph) [18]. 

The AO relation is satisfied for the limit of a Brownian random 
walk with a = 1 but shows a large deviation for fractional 

Brownian walks on GSC. d w characterizes the dynamic process of 

diffusive motion and d f describes the static geometrical scaling of 

the mass of the molecule. Thus, the ratio dJd w is a measure of 

molecular crowding. 



GFP-tagged RPA subunit in living HeLa cell by fixing 
a = 1 [14]. The corresponding average apparent 
translational diffusion coefficient for the fast mobile fraction 
was D = 0.153- 10~ 10 [m 2 /s]. The measurements were 

app L J 

performed with one-photon FCS [14]. In Fig. 6, we strictly 
observe compliance with the scaling laws of the simulated 

0.153- 10" 10 [m 2 /s] for 



apparent diffusion coefficient D 



all simulated a values over at least 3 orders of magnitudes of 
time scale T. The very small T values are below the 
experimentally accessible time scale because of the 
diffraction limit of the confocal optics used. As shown in 
Fig. (6A, C, E, G, I), the linear plots of simulated MSD(r) 
versus x clearly reveal translational anomalous motion for 
OC ^ 1 . For better comparison, normal diffusive motion is 
indicated by dashed lines. In Fig. (6B, D, F, H, J), the 
apparent diffusion coefficients D app (T) are plotted as lines 

from the double-logarithmic plot of the simulated MSD(r) 
versus X by means of Eqn. (3). Those lines give the scaling 
law of the measured D for different a values. Since in 

real experiments there are only a limited number of data 
points D (z) available, we verify the approximation for 
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which each single value of D app is found at time X from the 

slope of the straight line between the origin and MSD(x) [9]; 
this slope is 6 • D app (x) and the approximated D app values 

are given as dots in Fig. (6B, D, F, H, J). We obtain good 
agreement between the approximated D app values (dots) and 

the exact theoretical scaling laws (solid lines) in Fig. (6B, D, 
F,H,J). 

In Table 3, simulated data are given to account for the 
diffusion time T D 2 (x), e.g. of the slow mobile fraction of 

GFP-tagged RPA subunit in living HeLa cells [14], at 
different time scales X. The average diffusion time T D2 = 

5.340 ms over all time scales X was measured across the 3D 
observation volume AV of 0.14 fL by fixing a = 1 [14]. The 
measured value matches the expectation. Table 3 demon- 
strates that the diffusion time decreases for smaller values of 
a ^ 1 but increases for a larger time scale T at a given value 
of a ^ 1 . The effect of translational anomalous motion is 
substantially greater if a differs much from 1 . An a value 
close to 1 contributes little to the time dependence of 
subdiffusive motion. Thus, it is important to measure the 
anomalous dynamics on different time scales Tand to couple 
the analysis of how experimental parameters change with 
predictions from mechanistic models. 

So far we considered only one diffusing component. 
Now, we shall examine to what degree of accuracy 
translational anomalous motion can be distinguished in two- 
component systems. Usually, the measured diffusion time 
T D , the photon count rate detected from each kind of 

molecules, the relative quantum yield difference, and the 
molar concentration of each component in dilute solution are 
found by fitting an autocorrelation function for two different 
components to the experimental data [31]. In a situation of 
two labeled components in dilute solution, e.g. fast and slow 
mobile fractions, the measured diffusion time should differ at 
least by a factor of 2. A difference by a factor of 2 in 
diffusion times corresponds, for ideal spherical molecules, to 
a mass ratio of 1:8. Under certain conditions in dilute 
solution, a ratio of 1:5, e.g. for DNA/protein interactions, 
may be sufficient. In well-defined systems of dilute solution, 
differences in diffusion time down to 1:1.6 can be detected 
[31]. For globular molecules freely diffusing in dilute 
solution, this corresponds to a mass ratio of labeled ligand to 
formed complex or, generally speaking, of fast and slow 
mobile fractions of 1:4. Unfortunately, in a crowded 
environment of a living cell or its nucleus the motion of 
molecules are complicated as shown in Table 3 and the 
correlation data will be difficult to interpret. We present a 
criterion that rules out differences in diffusion times of two 
components which cannot be resolved in the correlation 
analysis of fluorescence fluctuation time series (photon 
streams) by FCS or FCCS or time-dependent fluorescence 
imaging under anomalous, subdiffusive motion. Assume that 
there is a fixed ratio of two diffusion times T D , ^ T D2 



corresponding to the masses of two components with m l ^ 
m 2 ■ If the resolution limit of T D 2 to T D , in the auto- or 
crosscorrelation function is given by, for example, 
N = 4 1/3 1.6 for comparable quantum yields of the two 
components, high fluorescence signal (number of photons 
detected for each species) and a molar fraction of the formed 
slower species that should not fall below 0.5 as has been 
shown by experiments and simulations under normal 
diffusive motion in dilute solution [31], the following holds 

^>K. (15) 

^ D,\ 

Our Eqn. (15) guarantees that T D , < X D 2 . In the case of 

two-color crosscorrelation or/and two-color fluorescence 
imaging (co-localization experiments), the second 
component is also a two-color species gr, for example a 
different molecular complex gr with an additional ligand. If 
we use this relation and additionally assume that 

T D 2 = T D j + At o with At^, the deviation from X D x 

satisfying our previous assumption of the resolution limit 
Eqn. (15), we then first derive a relation, which contains 

T D2 only and the deviation At D corresponding to a 

resolution limit for two different molecule species (two 
components) 

T D2 — -At d <0. (16) 

Our Eqn. (16) is valid for At™ < At b < At™ x . Here 
AtJT =T A2 -(K-1)/K and At™ =t™ x -(X -1)/R , where 
is the maximal value of the diffusion times observed in 
the set. This criterion allows us to select those values of 
T D2 , which satisfy Eqn. (15) in connection with the 

independent variables of the time scale T and the anomalous 
exponent (X (see Fig. 7, left panels). If we now change At 
we observe that only a fraction of the available values for 
T D 2 will satisfy relation (15) for At in a certain range. The 

selection of such resolution range in T D2 according to the 

mass ratio X = 4 1 3 — 1.6 reported for normal diffusive 
motion in dilute solution under well-defined and optimized 
experimental (optical) setup conditions [31] is shown in Fig. 
7 (right panels). These simulations demonstrate that the 
ability to distinguish between two different molecular 
species in a crowded environment depends strongly on the 
time scale T at which the subdiffusive motion occurs and on 
the degree of crowding quantitatively expressed by the OC 
value. Under anomalous subdiffusive motion due to 
molecular crowding, the specified value of the difference in 
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Fig. (6). Linear plots of simulated MSD(T) versus T (panels A, C, E, G, I) and corresponding double-logarithmic plots of simulated apparent 
diffusion coefficients -D app versus T (panels B, D, F, H, J) for different values of the simulated anomalous exponent a. The temporal 

autocorrelation defines the time scale Tof diffusion. In FCS, the autocorrelation is the crosscorrelation of the fluorescence signal with itself 
and is obtained by comparing a measured value at a time t with that at a time delayed by T. T is given in time steps covering a range from 
about 10 (J. s to 100 ms. For all simulated a values, we strictly observe compliance with the scaling law of the simulated apparent diffusion 
coefficient -D app = 0.153 • 10~ 10 [m 2 /s] of the fast mobile GFP-tagged RPA subunit in living HeLa cell over at least 3 orders of magnitudes 

of the time scale T. The very small rvalues are below the experimentally accessible time scale because of the diffraction limit of the confocal 
optics used. 
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Table 3. Simulation of the Time Scale Behavior of the Slow 
Mobile Fraction of GFP-Tagged RPA Subunit in 
Living HeLa Cells. The Average Diffusion Time 

Td,2(t) = 5.340 ms Over All Time Scales t Across the 
3D Observation Volume AV of 0.14 fL was Measured 
by Fixing a = 1 [14]. The Subdiffusive Motion in a 
Crowded Environment is Quantitatively Character- 
ized by the Anomalous Exponent a. The Simulated 
Time Scale T Covers a Time Range from about 10 (is 
to 100 ms at Which Diffusive Motion Occurred [14] 
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1.000 
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diffusion times of two components At b by a factor of 1 .6 at 

their mass ratio of 1:4 measured for normal diffusive motion 
under well-defined and optimized experimental conditions at 
comparable quantum yields for both components, a molar 
fraction of 0.5 and a high fluorescence signal [31] does not 
anymore guarantee that the two components can be resolved 
at all time scales T of the measurement and at all a values for 

In Fig. (7), we presented our simulated data for the 
resolution limit of two fractions with different molecular 
masses. We have first demonstrated that the crowding 
exponent a determines the resolution of differences in 
diffusion times in addition to the photophyscial parameters 
well-known for normal motion in dilute solution. For 
optimization of quantum yields of two components, number 
of photons detected for each species and the molar fraction 
of the formed slower species that should not fall below 0.5, 
etc., we refer to the experimental and simulated data 
available for normal diffusion in dilute solution [31]. Under 



optimized experimental conditions of the photophyscial 
parameters for two fluorescent components, we show that no 
better resolution of measured diffusion times can, in 
principle, be obtained than those given in Fig. (7) for a 
chosen K = 4 m — 1.6. This is due to the crowding of the 
intracellular environment. 

The argonaute (Ago) family of proteins are key 
regulators of RNA interference (RNAi) and RNA activation 
(RNAa) that function, in part, by recruiting small duplex 
RNAs [32-39]. It is well known that duplex RNAs [termed 
small interfering RNAs (siRNAs)] guide Ago proteins to 
complementary transcripts in the cytoplasm of cells to 
trigger RNAi [36-39]. Duplex RNAs [termed small 
activating RNAs (saRNAs)] can also activate gene 
transcription in the nucleus of cells to facilitate RNAa [32- 
35]. However, the mechanism of RNAa is complex and 
evidence is limited demonstrating physical interactions of 
Ago proteins with targeted nuclear sequences. We first apply 
our theoretical predictions under conditions of unbroken 
ergodicity to the time resolutions of molecular complexes 
formed in human prostate cancer cells (PC-3) transfected 
with fusions of green fluorescence protein (GFP) and 
argonaute proteins (Agol, Ago2) to determine the lower 
limit of time resolution according to mass differences. 
Fluorescence images show cellular distribution of GFP-Agol 
and GFP-Ago2 following overexpression in PC-3 cells (Fig 
8A). Note GFP-Agol and GFP-Ago2 can be found in both 
nuclear and cytoplasmic compartments of cells. GFP-fusion 
protein GFP-H2B was also transfected into PC-3 cells as a 
control for nuclear-specific localization; H2B is a core 
histone protein that functions exclusively in the nucleus (Fig 
8A). We believe that diffusive motions of any GFP-Ago 
protein (free GFP-Ago) bound at targeted sequence (GFP- 
Ago-Complex, bound GFP-Ago) can be time-resolved in 
living cells and their compartments if they fulfill the 
theoretical resolution limit (Fig. 8B). The molecular weights 
of both Agol and Ago2 are about 100 kDa, while the 
molecular weight of GFP is 27 kDa. The measured apparent 
diffusion coefficient of GFP in living cells is around 
0.275 • 10~ 10 [m 2 /s}. Based on these experimental values, we 
simulate the diffusion times and lower limits for the time 
resolution of two components in the fluorescence images. In 
the examples so far chosen, good agreement between the 
measured apparent diffusion coefficients at a = 1 and the 
known molecular weights show that our predictions are 
reasonable. In Fig. (8B), we present the simulated diffusion 

times T D2 for a GFP-Ago located at the target sequence 

(second component, bound GFP-Ago that is also called GFP- 
Ago-Complex) in the presence of 'free' GFP-Ago, i.e. 
unbound GFP-Ago (first component), for different simulated 
a values. Differences in diffusion times between the two 
components in a crowded environment can only be detected 
down to the X value for the chosen mass relation between 
GFP-Ago located at the target sequence (second component, 
bound GFP-Ago) in the presence of free GFP-Ago (first com- 
ponent) with K=4 1/3 = 1.6, K = 6 1/3 = 1.8, N=8 1/3 == 



538 Current Pharmaceutical Biotechnology, 2010, Vol. 11, No. 5 



Baumann et al. 




Fig. (7). The crowding exponent a determines the resolution of differences in diffusion times in addition to photophyscial parameters well- 
known for normal motion in dilute solution. Panel A: Total set of simulated diffusion times as function of X and a. The graph contains the 
behavior of diffusion times for some of the generators. It is obvious that x D decays with decreasing a and increases due to the scaling law 

T a in the direction of T. Panel B: Subset of simulated diffusion times 1 D satisfying the relation T fl2 — X'T fl |. Here, we 

choose K = 4 I/3 — 1.6. Panel C: Total set of diffusion times T D2 as function of Tand a. Simulated a values from top to bottom: 1.0, 0.977, 

0.924, 0.884, 0.829, 0.800, 0.745, 0.640. Panel D: This subset of values T D2 = T D l + At shows that for anomalous diffusive motion caused 

by molecular crowding the resolution limit given by the mass ratio t D2 /T Dl > K is not anymore guaranteed over all time scales T. 

Simulated CC values satisfying the chosen mass relation K = 4 - 1 .6 from top to bottom: 0.829, 0.800, 0.745, 0.640. 



2, X = 10 1/3 - 2.15, X = 100 1 ' 3 - 4.64. It is apparent from 
Fig. (8B) that the diffusion times of the second compound, 
Td,2 = Td,j + At ; At = Td,j • ( X -1), which can just be resolved 
in fluctuation time traces, decreases with decreasing a. This 
is predicted by the theory. In stationary FCS, the time spent 
dwelling on a spot between t = 0 and T D is partially wasted, 

since little new information is obtained in time intervals 
much smaller than the sample correlation time apart from 
overcoming photon statistics noise. This spare time can be 
effectively utilized by scanning into and observing new 
volumes, each of which are visited only occasionally. Single- 



molecule studies do not always use information from one 
individual molecule only or interpret every single molecule 
in the bulk phase [40]. Bridging single-molecule approaches 
with ensemble averages should yield interesting results [40, 
41]. The results assume that uncorrelated photons are 
measured which means they are statistically independent. 
Because the experimental apparatus requires time intervals 
of several milliseconds up to seconds and even longer, the 
photon correlations will be lost for times much longer than 
the coherence time that is the inverse of the bandwidth of the 
laser. This phenomenon can be understood by noting that if 
the simulation/measurement time T is very large, many 
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Fig. (8). Fluorescence images of PC-3 cells transfected with GFP-Agol or GFP-Ago2 overexpression vector. A: The images were provided 
courtesy of Dr. Long-Cheng Li at the University of California, San Francisco. GFP-H2B served as a control for nuclear-specific localization. 
Although enrichment varied, GFP-Agol and GFP-Ago2 were visible in the cytoplasm and nucleus of prostate cancer cells. B: We 
demonstrate that GFP-Ago complexes located at target sites must fulfill the resolution criterion N given in Eqn. (15). Simulated diffusion 
times T D2 f° r GFP-Ago located at the target sequence (second component, GFP-Ago-Complex) in the presence of GFP-Ago (first 

component) of living cells for different a values and chosen mass relations K = 4 1/3 — 1.6, X = 6 1 ' 3 — 1.8, K = 8 1 ' 3 — 2, K = 10 1 3 — 

2.15, K = 1 00 13 — 4.64 between putative GFP-Ago-Complex (second component) and GFP-Ago (first component). Simulated anomalous 
exponent OC , from top to bottom: 1.0,0.977,0.924,0.884,0.829,0.800. Measured values for o = 1 are in good agreement with known data of 
molecular weights of GFP-Ago2 and GFP-Ago 1 . The data points are from our theoretical predictions. 



fluctuations take place, and hence we measure an average 
value of the fluctuations by imaging and spectroscopy and 
not the fluctuation itself. The longer the time interval T, the 
closer the measured value approaches the mean value. As a 
consequence, the measured statistics approaches the 
uncorrected Poisson distribution. We provided a direct test 
of single-molecule trajectories in solution-phase by means of 
fluorescence fluctuation microscopy [10]. Our analysis 
moves beyond unphysical assumptions of diffusive 
measurements by time-dependent fluorescence fluctuation 
microscopy. 

The currently well accepted approach to measure a single 
molecule as it flows through a well-defined probe/ 
observation volume is often not true single molecule [42]. 
Although there is only one analyte molecule in the 
observation volume during the measurement, poor signal-to- 
noise requires that bursts from many analyte species must be 
averaged in order to achieve a reasonable signal-to-noise 
ratio. This makes it difficult to distinguish between rare 
confirmers with a strong signal that occasionally pass 



through the observation volume (or confirmers in dynamic 
equilibrium) from a mixture of stable confirmers (Richard. 
A. Keller, personal communication, Los Alamos); molecules 
that travel fast get there first and have less time to diffuse 
and the diffusion width is small. In very dilute solutions 
without flow, with very high probability the first molecule to 
enter the observation volume is the molecule that just left. 
The reentry time depends on the size of the observation 
volume, the diffusion coefficient and the molar bulk 
concentration of other molecules of the same kind that are 
not the original molecule [43]. Fluorescence molecule 
counting is, strictly speaking, a method to determine the 
experimental single-molecule regime for measuring just one 
individual molecule at a time [10]. It provides the theoretical 
and experimental criteria to verify not only the average 
molecule number but rather the molecular number of 
independently distributed molecules in dilute solution [10] as 
well as in crowded and highly heterogeneous living cells and 
their compartments (Baumann, Foldes-Papp, 2010, in 
preparation) [44]. 
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3.2. Translational Anomalous Motion with Broken 
Ergodicity 

Single biomacro molecules in in vitro systems measured 
by immobilization on cover slips display a heterogeneous 
molecular behavior. These are studies on one and the same 
single molecule, i. e. the selfsame molecule, over the whole 
measurement time T [42, 43]. However, we consider non- 
immobilized molecules [10, 43]. We present here the 
simulation results of translational subdiffusive motion under 
crowded condition where the time step becomes now a Levy 
distribution of waiting times [8]. To fulfill the experimental 
single-molecule regime the molecule number N per 3D 
observation volume AV has to be adjusted to N « 1 [10, 
43], and the single-molecule detection regime comes close to 
the average N = 0.048 and N = 0.0057 molecules and 
particles, respectively, per AV [45]. However, this 
experimental condition does not guarantee that one and the 
same molecule, i.e. the selfsame molecule, is measured over 
the whole measurement time T [10, 43, 46, 47]. A 
mathematical approach that rigorously solves the problem 
and produces the meaningful time for studying just one 
single molecule without immobilization or hydrodynamic/ 
electrokinetic focusing in solution and living cells is the 
meaningful-time concept first given in ref. [43, 46, 47]. 
These exact physical equations capture the reentry time of 
the selfsame molecule into the observation volume AV and 
depend on the size of the observation volume AV, the 
diffusion coefficient and the molar bulk concentration of 
other molecules of the same kind that are not the original 
molecule [43, 46, 47]. For the present results, we assume that 
the selfsame molecule is observed. Fig. (9) shows three 
examples for simulated data of the selfsame molecule from 
which MSD and the power law behavior can be evaluated. 
Random walks are usually carried out on regular lattices, 
which have the disadvantage that the spatial and temporal 
steps are correlated. However, randomness in cellular 
systems can be either spatial (anomalous) or temporal 
(heterogeneous). In order to separate both processes we 
introduced in the previous sections anomalous random walks 
on fractals that represented crowded environments. In order 
to mimic heterogeneous behavior in crowded environments 
of living cells, we need to introduce so-called continuous 
time random walks (CTRW). CTRWs were originally 
performed on regular lattices [48]. In this section, we 
combine both ideas and perform random walks on fractals 
using CTRWs [49]. This purely stochastic molecule behavior 
leads to subdiffusive motion with broken ergodicity in our 
simulations. 

In a CTRW, the molecule has to wait for a time t on each 
site of the fractal before performing the next step. This 
waiting time is a random variable independently chosen at 
each new step according to a continuous distribution i^(-r) ■ 
In our case, a stable Levy distribution with Levy index 3/2 is 
used for the waiting time steps ^( T ) = e -"< 2l > /(V^x 3 ' 2 ) • The 

first and second moments of this distribution do not exist. 
We assume that the molecule can step to one of the nearest 
neighbors generated by the fractal generator. In addition, 
y/(t) is independent of the location of the molecule, i.e. we 

decouple spatial and temporal coordinates. One can think of 



this process as a diffusive motion among traps, but the 
trapping time will change for each site on the lattice if this 
site is visited again during the random walk. If we perform a 
walk on a regular lattice, disorder can be mimicked through 
waiting time distributions y/(z) displaying long-time tails 

y/-(T) ~ T" <1+7) with 0<7<1. (17) 

The probability to perform n steps during time T is 
denoted by % n (t) , which is related to the waiting time 
distribution by the Laplace transform 

L{Xn(*)) = Zn(s) = (VW • C 1 - V(*)) / s. (18) 

This probability is needed to analyze the MSD for a 
random walk on a fractal support carried out as a CTRW [49, 
50]. The MSD traveled by the molecule during the time T is 
given by 

(?(T))=X^yZn(*), (19) 

/- 2 \ 

where lr„) is the average distance traveled in n steps on 

/- 2 \ 

the fractal. For fractals, we found that lr„)~n a so that 

after carrying out the summation in the continuous limit the 

/- 2 \ 

asymptotic behavior ( r „ ) is found to be 

(7(T)j~T ay =T\ (20) 

We observe in our simulations that, contrary to what is 
expected from random walk on fractal support (RWF), the 
temporal averaged MSD r Eqn. (10) leads to a simple 
diffusive behavior with scaling exponents that strongly differ 
from one track to another. This distribution of scaling 
exponents renders a system inhomogeneous: an ensemble of 
simple diffusers with different diffusion coefficients. Since a 

and f satisfy 0 < a,f < 1 , so does the product a ■ y and the 
diffusive motion is thus subdiffusive. The main point here is 
that the ensemble average and the time average do not 
deliver the same scaling exponent anymore, i.e. we have 
broken ergodicity. 

Since the experimental conditions to measure a selfsame 
molecule over an extended period of time, at which biology 
is taken place, in living cells or even in dilute solution-phase 
are very restrictive [43, 46, 47], we need to perform the time 
average over a subpopulation of different single molecules of 
the same kind: 




Our experimental single-molecule regime given by Eqn. 
(21) differs from averaging over the whole molecule 
ensemble suggested by Meroz, Sokolov, Klafter (2010) [8]. 
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Fig. (9). MSD for a specific generator of the fractal grid (top left). Top right the apparent diffusion coefficient is shown for the random walk 
on fractal support (RWF). For each track 10 4 steps on the grid were performed in a random walk on fractal support and in a CTRW on the 
fractal support. The ensemble for the RWF consist of 10 4 tracks from which the MSD follows (bottom left). For the CTRW the different 
tracks all have the same number of steps. The CTRW with the selected fractal generator shows a broad distribution of fractal dimensions. 
However, the ensemble mean for the same fractal generator shows a single reproducible fractal dimension of a = 0.830. The specific fractal 
dimensions following from the temporal averaging in the MSD are y = {0.751,0.821,0.935} shown in the lower right figure. Since the RWF 

scaling exponent is fixed to a = 0.830 we can derive the time scaling from the relation y = a-y to be y = {0.905,0.989,1.126} ■ This behavior 

shows that for the MSD the ensembles average and the time average is not any more the same; i.e. we have broken ergodicity. 



For time averages over subpopulations of single molecules 
of the same kind, the temporal auto- and crosscorrelation 
functions are found for the first time 

G(.) = l.(l + (W TD ff.(l + (l/^).(r/ Tfl fp 2V2 + DC, (22) 

y = ay, (23) 

Here, T D is the diffusion time and T D is the two-color 

cross-correlated diffusion time given by Eqn. (2). s is the so- 
called structural factor defined as s = (Oj(O x _ . In the case 

of 3D diffusive measurements, dim = 3. For 2D diffusive 
measurements (e.g., in membranes), dim = 2. DC is the 
limiting value of G(t) for t — > °° , which is normally 1. By 

fitting Eqn. (22) to the time traces that are recorded for the 
subpopulation of single molecules of the same kind without 
interacting partner, e.g. without ligand, in the crowded 
environment of a living cell and their cellular compartments, 
respectively, y = a is obtained. Knowing the crowding 

parameter a for the cell type and cellular compartment type, 
the heterogeneous parameter y can be extracted from the 
measurements in the presence of the interacting reaction 
partner, e.g. ligand, for the same a. y (Eqn. (23)) is not a 
simple fitting parameter in the temporal auto- and two-color 
crosscorrelation functions (Eqn. (22) and Eqn. (2)) since it is 



related via Eqn. (23) to the physical model given in Eqn. (19) 
for our proposed Eqn. (21). We already derived an analytical 
solution for y in the special case of y = 3/2 [44]. 

CONCLUSIONS 

The generally applied assumption of normal Brownian 
motion of macromolecules works well for in vitro 
experiments under dilute conditions, but does not apply for a 
crowded environment of living cells and their cellular 
compartment, e.g. the nucleus. Fixing a = 1 in auto- and 
crosscorrelation functions for one and two species (compo- 
nents) results only in a coarse average over all measured 
time scales X of subdiffusive motion. Our simulations show 
that determination of molecular weights from measured 
diffusion times and apparent diffusion coefficients, respec- 
tively, in temporal auto- and two-color crosscorrelation 
analyses as well as from fluorescence imaging data are 
difficult to interpret and biased in crowded environments of 
living cells and their cellular compartments unless the 
anomalous dynamics on different time scales T is coupled 
with the analysis of how experimental parameters change 
with predictions from simulated diffusion times and apparent 
diffusion coefficients, respectively, and mechanistic models. 

The basic property of ergodic systems is that the time 
average and the ensemble average on the tracks result in the 
same MSD. Nonergodicity is given if the time average and 
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the ensemble average on the tracks are different. The noner- 
godic situation is strictly valid only for the selfsame mole- 
cule. A mathematical approach that produces the meaningful 
time for studying the selfsame molecule (just one single 
molecule) without immobilization or hydrodynamic/ electro- 
kinetic focusing in solution and living cells is the meaning- 
ful-time concept first derived in ref. [43, 46, 47]. These exact 
physical equations capture the reentry time of the selfsame 
molecule in the observation volume AV and depend on the 
size of the observation volume AV, the diffusion coefficient 
and the molar bulk concentration of other molecules of the 
same kind that are not the original molecule [43, 46, 47]. 
Since the experimental conditions to measure a selfsame 
molecule over an extended period of time, at which biology 
is taken place, in living cells or even in solutions are very 
restrictive, we need to perform the time average over a 
subpopulation of different single molecules. The temporal 
auto- and crosscorrelation functions are found for the first 
time with the crowding parameter a for the cell type and 
cellular compartment type, respectively, and the heteroge- 
neous parameter y that can be obtained from the measure- 
ments in the presence of the interacting reaction partner, e.g. 
ligand, with the same a value. The product y of the 
crowding parameter a and the heterogeneous parameter y is 
not a simple fitting parameter in the temporal auto- and two- 
color crosscorrelation functions because it is related to the 
proper physical models of anomalous (spatial) and 
heterogeneous (temporal) randomness in cellular systems. 

This work extends previous work on trafficking of 
mature miRNA-122 into the nucleus of living liver cells, 
where we first presented FCS/FCCS-artifact free subdif- 
fusive measurements of the formed fluorescent complex in 
the cytoplasm and nucleus by two-photon imaging [51]. A 
similar biological interpretation of miRNA traveling into the 
nucleus was proposed independently by Ohrt et al. [13, 52] 
from quantitative two-color FCCS measurements and respec- 
tive calculations of molecular weights in living cells but 
without taking into account the difficulties and artifacts in 
the proper interpretation of quantitative subdiffusive mea- 
surements that we presented above. In crowded environment 
of living cells and their compartments, diffusive motion is 
characterized by the time-dependence of the mean-square 
displacement and the apparent diffusion coefficient; the 
measured diffusion time is not anymore a fixed constant at 
all time scales X for a given size of the observation volume. 
Here, we first show that plausible biological mechanisms 
from FCS/FCCS and fluorescence imaging data are highly 
questionable without proper quantitative physical models of 
subdiffusive motion and temporal randomness. At best, such 
quantitative FCS/FCCS and fluorescence imaging data are 
difficult to interpret in living cells under crowding and 
heterogeneous conditions. It is challenging to translate pro- 
per physical models of anomalous (spatial) and heteroge- 
neous (temporal) randomness in living cells and their cellular 
compartments like the nucleus into biological models of the 
cell biological process under study testable by single- 
molecule approaches in living cells. Otherwise, quantitative 
FCS/FCCS and time-dependent fluorescence imaging mea- 
surements in living cells are not well described and cannot 
be interpreted in a meaningful way. 
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